Structure and magnetic properties of nanocrystalline PrCo 3 
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The structure and magnetic properties of nanocrystalline PrCo3 prepared by high-energy milling 
technique have been investigated by means of X-ray diffraction using the Rietveld method coupled to 
Curie temperature and magnetic measurements. The as-milled samples were subsequently annealed 
in temperature range from 750 to 1050° C for 30min to optimize the extrinsic properties. From x-ray 
studies of magnetic aligned samples, the magnetic anisotropy of this compounds is found uniaxial. 
The Curie temperature is 349 K and no saturation reached at room temperature for applied field of 
90kOe. The coercive field of 55kOe and 12kOe measured at 10 and 293 K respectively is obtained 
after annealing at 750° C for 30 min suggests that nanocrystalline PrCo3 are interesting candidates 
in the field of permanent magnets. We have completed this experimental study by simulations in 
the micromagnetic framework in order to get a qualitative picture of the microstructure effect on the 
macroscopic magnetization curve. From this simple model calculation, we can suggest that the after 
annealing the system behaves as magnetically hard crystallites embedded in a weakly magnetized 
amorphous matrix. 

PACS numbers: 75.50.Bb, 75.50.Tt, 76.80.+y 



I. INTRODUCTION 

The intermetallic compounds, RC03 (R= rare earth 
element) have attracted great interest caused mainly by 
their wide applications as permanent magnets. These 
compounds present excellent magnetic properties such as 
large magnetocrystalline anisotropy and important satu- 
ration magnetization at room temperature. This inter- 
esting magnetic performances of the intermetallic com- 
pounds formed by alloying rare earth and 3d transition 
metals are due to the combination of the complemen- 
tary characteristics of 3(i-itinerant and 4/-localized mag- 
netism. The coercivity is one of the most important ex- 
trinsic properties of a permanent magnet. This prop- 
erties +depends on the microstructure of the material. 
Recently nanocrystalline compounds have been exten- 
sively investigated, their aim is to find the appropriate 
nanocrystalline state for a highest coercive field. 

The most employed process to prepare these materials 
are rapid quenching by melt spinning [l[ and mechanical 
alloying Q. More recently high-energy milling technique, 
has also been used to produce high coercivity nanocrys- 
talline powders. Rivera-Gomez et al are reported in Q 
a high coercivity for PrCos compounds milling and an- 
nealed at 800°C for 1 min with H c = 17.3 kOc. A high co- 
ercivity of Sni2Coi7 nanocrystalline compounds has been 
obtained by using this method [J] where H c = 9.6 kOe 
after milling and annealed at 800° C for 30 min. 

In this paper, we report the elaboration of nanocrys- 
talline PrCo3 by high-energy milling, and their structural 
and magnetic properties in amorphous state and after re- 
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crystallization. Hard magnetic properties of the as milled 
PrCo3 have been improved by a controlled nanocrystal- 
lization, which optimize the size of grains in order to 
obtain a high coercivity. 

We also propose a simple model for the system after 
annealing based on a simple lattice of crystallites embed- 
ded in an amorphous matrix. The simplest realization of 
such a model is considered here, in order to bring out at 
the qualitative level the microstructure effect on the mag- 
netization curve. This model is treated in the framework 
of micromagnctism, i.e. through a continuous medium 
type of approach, and the calculations are performed by 
using the MAGPAR @ code based on the finite elements 
method. 



II. EXPERIMENT 

A polycrystalline PrCo3 was prepared by melting high- 
purity starting elements (> 99.9%) in an induction fur- 
nace under a protective environment (Ar atmosphere). 
Their homogeneity was insured with five times consec- 
utive melting. These ingots were used as pre-alloys to 
manufacture samples by high-energy milling. All pow- 
der handling was performed in a glove box under a high- 
purity Ar atmosphere with an O2 and H2O rate around 
1 ppm with vial sealed glove box. 

The mixtures were scaled in hardened steel vial to- 
gether with 15 mm diameter steel balls. Milling was per- 
formed for 5 h under a high purity Ar atmosphere with a 
ball to powder ratio of 15 : 1, in a high-energy planetary 
ball mill Fritsch P7. These milling conditions correspond 
to kinetic shock energy, shock frequency and injected 
shock power values respectively equal to 0.81 J/hit, 62 Hz 
and 19.5watt/g. After milling, the powder mixtures was 
wrapped in tantalum foil and sealed in silica tubes under 
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a vacuum of 10~ 6 Torr, then annealed for 30 min at tem- 
perature between 750 and 1050°C followed by quenching 
in water. 

The purity of the final phase was checked by X-ray 
diffraction (XRD) using Bruckcr diffractomctcr mounted 
with the Bragg-Brcntano geometry and automatic diver- 
gence slit, using the CuKa radiation. The unit param- 
eters were measured with Si as standard (a— 5.4308 A) 
to insure a unit-cell parameter within the accuracy of 
±lxl0 -3 A. The counting rate was 22 s per scanning step 
and step size of 0.04°. To determine the microstructural 
alloy properties (crystallite size) and crystallographic 
properties (phase identification and lattice parameters), 
XRD patterns have been analyzed by Rietveld method. 

The data treatment was carried out by a Rietveld re- 
finement @,J3 as implemented in the FULLPROF com- 
puter code [a with the assumption of a peak line pro- 
file of Thompson-Cox-Hastings allowing multiple phase 
refinement of each of the coexisting phases • It takes 
into account the broadening of the diffraction lines and 
gives, via the Scherrer formula, the autocoherent diffrac- 
tion domain size. The weight percentage of each of the 
coexisting phases is directly obtained with the approxi- 
mation of a Brindlcy absorption factor equal to 1 . The 
goodness-of-fit indicators Rb and x 2 are calculated as 
usual and described earlier [n], [l2| ■ 

The magnetic anisotropy of PrCo3 was studied us- 
ing X-ray diffraction. The X-ray diffraction pattern was 
recorded on powder sample (grain size < 50 /im) oriented 
under an external magnetic field. The magnetic field was 
oriented perpendicular to the surface of the sample and 
the grains are oriented and shape the needles parallel to 
the field applied, the powder so directed is analyzed by 
XRD at room temperature, thus the x-ray pattern gave 
information about the easy magnetization direction. 

Curie temperature T c was measured on a differen- 
tial sample magnetometer MANICS in field of 0.1 kOe 
and heating rate of lOK/min with around 10 mg sample 
sealed under vacuum in silica tube in order to prevent 
oxidation under heating. The Curie temperature was de- 
termined from the M — T curves by extrapolating linear 
part of the M — T curves of the intersection with the ex- 
tended base line. The M — H curves (Magnetic hysteresis 
measurement) were plotted at T= 10 and 293 K with a 
Physical Properties Measurement System (PPMS) Quan- 
tum Design equipment and a maximum applied field of 
90 kOe on sample in epoxy resin. 



III. RESULTS AND DISCUSSION 

A. Structure analysis 

The PrCo3 compounds crystallize in rhombohedral 
PuNi3-type structure (space group: R3m). The crys- 
tal lattice of PrCo3, can be regarded, approximately, as 
consisting of the consecutive arrangement of alternating 
blocks of PrCos and PrCo2 along the c axis with the ra- 



tion 1 : 2 (Fig. Q}, having almost the same structure as 
in the corresponding compounds PrCos and PrCo2- The 
unit cell contains two non-equivalent crystallographic site 
for Pr ions, 3a and 6c, and three sites for Co: 36, 6c and 
18h. 

The X-ray diffraction (XRD) of the PrCo3 shows 
that the compound is single phase with the rhombohe- 
dral PuNi3-type structure (space group: i?3m). Fig. [5] 
presents the analysis result of XRD pattern of PrCo3 in- 
got using the Rietveld analysis. The lattice parameters 
of PrCo3 compound are a= 5.069 A and c= 24.795 A. 
The atomic positions, Rb, X 2 factors from Rietveld fit 
are given in Table 1 . These values of structural parame- 
ters are in agreement with the results obtained previously 
by [Hj]. This ingot is an appropriate precursor for us- 
ing as pre-alloys to manufacture samples by high-energy 
milling. 




FIG. 1: Crystal structure of PrCo3. Pr atoms occupy the 3a 
and 6c sites, and Co atoms occupy the 36, 6c and 18/i sites. 
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FIG. 2: Rietveld analysis for PrCo3 ingot. 
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TABLE I: a and c cell parameters, atomic positions, Kb and 
X 2 factors from Rietveld fit. 

a(A) 5.0690 
c(A) 24.7950 
z(6e)Pr 0.1414 
z(6c)Co 0.3337 
x(18/i)Co 0.5002 
y(18/i)Co 0.4998 
z(18/i)Co 0.0829 
c/a 4.8915 
V(A) 3 551.42 
Rb 5.050 
X 2 2.300 



Figure [3] shows the x-ray diffraction pattern of three 
samples of PrCo3, the as- milled powder and annealed at 
750 and 1000°C for 30min. 

The as-milled powder is typical of a quasiamorphous 
phase and does not show characteristic lines of diffrac- 
tion of the crystallized state, however, we observe in the 
low angles a flattened and large line around 42°. Such 
quasiamorphous matrix suggests the presence of various 
types of small crystallites. The diffraction peaks become 
more pronounced with increasing annealing temperature. 
The annealed samples at 750 and 1000°C for 30min re- 
veal the presence of PrCo3 phases crystallized. We have 
observed a small amount (around 0.5 wt%) of Pr203 in 
all samples annealed after milling. 

The mean diffraction crystallite size of the powder af- 
ter annealing at 750 and 1000°C for 30 min obtained by 
fitting the XRD diagram with Rietveld method, ranges 
between 31 and 45 nm respectively, this implies that syn- 
thesized materials are noncrystalline. 
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FIG. 3: XRD pattern of the PrCo3 as-milled and annealed 
samples at 750 and 1000° C for 30 min. 



B. Intrinsic magnetic properties 

The Curie temperature, also called the magnetic order- 
ing temperature is a direct measure of the exchange in- 
teraction, which is the origin of ferromagnctism. This in- 
teraction depends markedly on the interatomic distance. 

Generally, the Curie temperature in rare-earth 
transition-metal intermetallic compounds is governed by 
three kinds of exchange interactions, namely, the 3c? — 3c? 
exchange interaction (JcoCo) between the magnetic mo- 
ments of the Co sublatticc, 4/ — 4/ exchange interac- 
tion (Jp r p r ) between the magnetic moments within the 
Pr sublattice, and the intersublattice 3d — 4/ exchange 
interaction (Jp r co)- The interactions between the rare- 
earth spins 4/ — 4/ are assumed to be weak and negligible 
in comparison with the other two types of interactions 
[p4| . There are two exchange interactions in PrCo3 com- 
pounds, positive and negative. When the distance of the 
Co-Co pairs is smaller than 2.45 A, the exchange interac- 
tions is negative, whereas at larger Co-Co distances the 
interaction is positive. 

As an example, the thcrmomagnetic magnetization 
curve of PrCo3 high-energy milling sample annealed at 
750°C for 30 min is presented in Fig. 2) This sample has 
a unique magnetic phase transition at the Curie temper- 
ature of 349 K. This result is in agreement with the data 
given by Lemaire [l5| . 

The low Curie temperature in PrCo3 compounds is 
due to the short Co-Co interatomic distances at the sites 
6c — 18h and at 18 h — 18h (Table 2), where the Co atoms 
couple antiferromagnetically. In this structure, the 6c — 
18h and at 18h — 18h interactions are strongly negative. 
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FIG. 4: Thermomagnetization curves of PrCo3 obtained after 
high-energy milling and annealed at 750° C for 30 min. 

Isothermal magnetization curves M — H obtained at 
293 K is represented in Fig. [5] for PrCo3. This figure 
shows that at 293 K the saturation is not reached, which 
gives evidence for a magnctocrystallinc anisotropy among 
the highest known in those kinds of systems. 

The pattern obtained on sample oriented under an 
external magnetic field, applied perpendicularly to the 
plane of the sample, shows only (00Z) Bragg peaks 
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TABLE II: Interatomic distance (A) for PrCo 3 



Site 


Co (36) 


Co (6c) 


Co (18ft) 


Pr (6a) 


Pr (6c) 


Co (36) 


5.07 


4.12 


2.54 


2.99 


5.07 


Co (6c) 


4.12 


2.92 


2.42 


4.56 


2.92 


Co (18ft) 


2.53 


2.42 


2.41 


2.92 


3.50 


Pr (6a) 


2.99 


4.56 


2.99 


3.18 


3.50 


Pr (6c) 


2.99 


2.92 


3.26 


3.50 


5.07 
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FIG. 5: Magnetization curves of PrCo3 obtained after high- 
energy milling and annealed at 750° C for 30min, measured 
at 293 K. 
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FIG. 7: Hysteresis loop of PrCo3 milling before and after 
annealing at 750° C for 30min, measured at 293 and 10 K. 



(Fig. [BJ), which amount to suggest that the anisotropy 
strengthens the peaks of diffraction of type (00/) and 
at room temperature the easy magnetization direction 
is parallel to the c-axis. A similar result is obtained for 
YCoz- x ¥e x by obtained by oriented powder under 
extern field. 




FIG. 6: XRD pattern of PrCo 3 oriented. 



C. Extrinsic magnetic properties 

Concerning the extrinsic magnetic properties, we have 
measured hysteresis loops at 10 and 293 K for all sam- 
ples annealed at different temperature (750 to 1050°C 



for 30min) and plotted M — H curves. Fig. [7] show 
as an example the hysteresis loops measured at 10 and 
293 K for sample annealed at 750°C. This sample ex- 
hibits high coercivities, we have obtained, for T = 10 K, 
H c = 55kOe, Mr = 56uem/g and M R /M max = 0.75. 
For T= 293 K, H c = 12kOe and Mr = 44uem/g and 
MR/Mmax = 0.60 inherent to nanocrystalline state. The 
coercivity of 55kOe at 10 K confirms the large magne- 
tocrystallinc anisotropy. 



IV. SIMULATION 

In order to get a first picture of the link between the 
structure at the nanometric scale and the extrinsic mag- 
netic properties at least at the qualitative level, we have 
performed micromagnetic simulations on a model for our 
system. Indeed, micromagnetic simulations provide a 
useful tool for this purpose and have been done for a 
wide range of nanocomposite magnetic systems [l7l - [22^ . 
It is worth mentioning that quite generally such simu- 
lations deal with models for heterogeneous systems in- 
cluding magnetic grains of different materials, typically 
hard and soft materials, the purpose being to character- 
ize the exchange |l7H20l| or the magnetostatic coupling 
between phases [23| ■ The general finding is that the inter- 
grains exchange coupling leads to a decrease of the coer- 
civity and an increase of remanence when compared to 
the values of the hard material component. The magnc- 
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tostatic coupling may be more complex due to the long 
ranged and the anisotropy nature of dipolar interactions 
pij . Here, we model a biphasic system, the phases of 
which are of the same chemical composition. Indeed, as 
we have seen in the experimental section, the as milled 
sample is amorphous, and then after annealing well crys- 
tallized grains whose size is controlled by the annealing 
temperature are formed. Our purpose, rather than to 
provide a systematic investigation, is to draw out some 
of the salient features concerning the sample structure 
which may explain the shape of the magnetization curve. 
We model the experimental sample as an assembly of 
well crystallized grains, represented by either spherical 
or cubic inclusions, embedded in the matrix made of the 
amorphous phase corresponding to the as milled mate- 
rial. The inclusion shape should be discriminated on 
physical grounds related to the crystallization growth at 
the temperature considered. This is beyond the scope 
of the present work, and we will be lead to discriminate 
the grain shape from the agreement with experimental 
results. We emphasize that the main difference between 
these two cases is the larger relative importance taken by 
the amorphous matrix, playing the role of grain bound- 
aries, in the model with spherical inclusions. 

Since the crystallites are well crystallized, they are 
characterized by the magnetic constants of the perfect 
massive crystal. We take the experimental value of the 
saturation magnetization per unit volume, J s = 0.723 T 
as deduced from both the saturation magnetization 
per unit formula [l5j |. namely and the volumic 

mass as obtained from the measured volume of the unit 
cell (see table fl} which corresponds to Z = 9. The 
value of the exchange constant is set equal to that of 
the archetype of hard material (Nd2Fci4B), namely 
A ex = 7.7 x 1CP 12 J/m and the magnetocrystallinc 
anisotropy constant is fitted on the experimental value 
of the coercive field we have got at T= 10 K, with the 
result K\ = 3.05 10 6 J/m 3 . The amorphous matrix is 
on the one hand non anisotropic [K\ = 0), and on 
the other hand characterized by a very small value 
of the exchange constant, set somewhat arbitrarily to 
A ex = 3.8 x 10 _13 J/m. Indeed, since the coercivity 
remains very high we expect in our two phase model only 
a negligible exchange coupling between the matrix and 
the inclusions. Then, concerning the matrix saturation 
magnetization, we expect a much smaller value than for 
the crystallized phase of the inclusions [25j]; therefore 
we consider either a non magnetic (J s = 0) or only 
weakly polarized (0.05 T < J s < 0.25 T) matrix. As we 
shall see these two last points can be justified by the 
magnetization curve shape. The inclusions have all the 
same size, while their easy axis, hi are distributed either 
randomly or in a solid angle < ©a/ centered on the 
direction of the external field, h ext in order to model the 
effect of the texturation by the field. The influence of the 
preferential orientation has been investigated only in the 
case of spherical inclusions. The inclusions are located 
on the sites of simple cubic lattice. The inclusion packing 



fraction is therefore limited to values smaller than ir/6 
in the case of spherical inclusions; we set the nearest 
neighbor distance d = 2.106 R close to the minimum 
value d = 2 R, which leads to <p = 0.445. We have not 
such limitation when dealing with cubic inclusions. In 
this later case the inclusions are separated by a grain 
boundary of thickness, say <5, and the volume fraction is 
then tp = (1 + 6/2R)~ 3 where the inclusion hedge length 
is set to 2 R. In the framework of the micromagnetism, 
the value of 6 cannot be reduced at will since 8 must 
not reach the microscopisc scale; here we consider either 
5/d = 0.10 or S/d — 0.237 where d is the minimum 
intergrain distance which leads to ip = 0.73 or 0.445 
respectively The embedding matrix is a parallepipcdic 
box including 4 layers of 64 lattice sites and the total 
number of particles taken into account is thus 256. The 
applied field is set normal to the layers. In order to 
simulate as far as possible the experimental case, we 
have chosen R= 17nm for the spherical inclusions and 
2i? = 32.5nm for the cubic inclusions. Accordingly 
the size of the matrix is 288 nm x 288 nm x 144 nm 
and in between 297 nm x297nm xl58nm and 
340 nm x340nm xl70nm in the case of spherical 
and cubic inclusions respectively The micromagnctic 
calculations are performed by using the MAGPAR code 
based on a finite elements scheme, from either the 
LLG equation or the total energy minimization. Most 
of the calculations have been performed with a mesh 
including 1632000 tctrahedra, and we have checked on 
a refined mesh of 3159000 tctrahedra the convergence 
versus the element size of the results in the case of 
spherical inclusions. 

The results for the model with spherical inclusions are 
displayed in figure ©. First of all, we consider the case 
of the non magnetic matrix (J s = 0) with randomly ori- 
ented particles. As expected, the coercivity takes a large 
value, namely totally governed by that of the particles, 
since the exchange constant of the matrix is too low for 
the particles to be exchanged coupled to and via the ma- 
trix. We have checked that even with J s = a significant 
value of the exchange constant for the amorphous matrix 
leads to a decrease of the coercivity, and this as already 
mentioned justifies our choice of a very small value of 
A ex for the matrix. However, the main qualitative char- 
acteristic which is missing in this first calculated magne- 
tization curve when compared to the experimental one is 
the small kink in the vicinity of zero external field. Such 
a kink may be the signature of an uncoupled two phase 
magnetic behavior. Hence, we have increased the value of 
the matrix saturation magnetization in order to get such 
a behavior and we find that J s = 0.05 T leads to a quali- 
tative agreement with the experimental curve. When the 
matrix J s value increases the kink amplitude becomes 
much larger than the experimental one. As a matter of 
fact, the kink relative amplitude behaves roughly as 

(1 - ip) jj mat) 

(1 - ip) + pdi inCl) 
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FIG. 8: Calculated de-magnetization curve in terms of the 
external field for the systems of spherical inclusions with 
tp = 0.45. Randomly distributed axes and matrix saturation 
magnetization J s = 0, dashed line; J s = 0.05 T, solid line. 
The whole hysteresis curve in this later case is obtained by us- 
ing the symmetric of the de- magnetization curve. J s = 0.05 T 
and Om = n/2.3, crosses; 8m = 7r/15, dotted line. 

The agreement between the calculated hysteresis curve 
and the experimental one, for randomly distributed axes 
and jj ma *' = 0.05 T is qualitatively satisfying. We 
have to mention however, that we display in figure 
© the magnetization curve reduced by the total sat- 
uration magnetization, M s , of the system the value 
of which, M s = ipJs inc ^ + (I — ip)Js mat ^ is significantly 
smaller than the experimental one. This is due both to 
the small value of the saturation magnetization of the 
amorphous matrix obtained to reproduce the amplitude 
of the kink at H ex = 0, and to the large fraction of the 
total volume occupied by the matrix. The model with cu- 
bic inclusions can remedy to this drawback since then the 
inclusion packing fraction is not limited as already men- 
tioned. The results obtained with tp = 0.73 are displayed 
on figure © for the matrix saturation magnetization 
J™ a *= 0.05 T, 0.125 T and 0.250 T. The first conclusion 
is that at the qualitative level, the hysteresis curve looks 
like that of the system with spherical inclusions which 
may be the consequence of the easy axis random orienta- 
tion. Now, since the matrix volume fraction is strongly 
reduced, the amplitude of the magnetisation kink due to 
the matrix magnetization reversal compares to the ex- 
perimental one with a larger value for J™ ai than in the 
model with spherical inclusions, namely J™ a * ~ 0.125 T. 
However, increasing the matrix saturation magnetization 
leads to an important matrix/inclusion magnetostatic 
coupling resulting in a non negligible reduction of the 
coercivity (see figure [9]). Thus the best agreement with 
the experimental results is obtained with the model with 
cubic inclusions, tp = 0.73 and jj" 10 *' = 0.125 T {i.e. 
j(mat) J j(mci) _ q^72). With this set of parameters, we 
get for the saturation magnetization of the composite sys- 
tem ji* ot) = 0.777 ji mc °. The simulated magnetization 
curve of the model with cubic inclusions and tp = 0.445 is 



quite close to the that of the model with spherical inclu- 
sions and is thus not shown. We also get an important 
reduction of the coercivity with the increase of the matrix 
saturation magnetization. 
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FIG. 9: Calculated de- magnetization curve in terms of the ex- 
ternal field for the systems of cubic inclusions with tp = 0.73. 
Randomly distributed axes and matrix saturation magneti- 
zation J 3 = 0.05 T, dashed line; J s = 0.125 T, solid line and 
J s = 0.25 T, dotted line. 

Then, and only for the model with spherical inclusions, 
we have examined the effect of the texturation via the 
easy axes distribution. In this calculation, the maximum 
value of the polar angle of the axes ri, is either 7r/15 or 
7r/2.3 for a strong or a very weak preferential orientation 
respectively. Notice the strong sensitivity to the varia- 
tion of 0m! indeed the random distribution corresponds 
to ©m = tt/2. We emphasize that this orientation effect 
may represent a consequence of a non spherical shape for 
the particles. Finally, we have also calculated the magne- 
tization curve for particles of size R = 35 nm with nearly 
the same result at least in the case on a non magnetic 
matrix. Therefore, as expected, we conclude that the 
strong size effect observed on the coercivity of nanocom- 
posite materials is mainly due to the dependence of the 
crystallite magnetic properties on the size. 



V. CONCLUSION 

In summary, the PrCo3 intermetallic compound, crys- 
tallize in rhombohcdral PuNi3-type structure (space 
group: R3m). This compound is a highly anisotropic uni- 
axial fcrromagnet with the easy magnetization direction 
is parallel to the c-axis. No saturation reached at room 
temperature for applied field of 90 kOc. High coercivi- 
ties of about 12kOe at 293 K and 55kOe at 10 K can be 
reached in PrCo3 alloy prepared by high-energy milling 
method followed by annealing at 750°C for 30min. High 
energy milled method following the controlled nanocrys- 
tallization brings high magnetic properties with smaller 
grain size more adequate to favor an available microstruc- 
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ture to perform this properties, so these results show that 
PrCo3Compounds are promising materials for permanent 
magnet applications. From the model calculations of the 
macroscopic extrinsic properties we can propose a simple 
picture of the microstructure. Indeed, the magnetization 
curve presents a signature of a system made of two un- 
coupled phases: the magnetically hard crystallites and 
the amorphous matrix. The presence of the small kink 
in the magnetization curve is analyzed as the reversal of 
the matrix magnetization; the localization of this kink at 
nearly zero applied field and the its small amplitude are 
indications that the matrix is, as expected, magnetically 
very soft and characterized by a small saturation magne- 
tization. This latter point must be taken with care, since 



here we propose a simple lattice for the crystallite which 
may lead to an overestimated fraction of the total volume 
occupied by the matrix also with the model with cubic 
inclusions. It is important to notice that when dealing 
with the reduced magnetization, M/M s we find that the 
influence of the grain shape is rather weak. 
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